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Abstract 

We propose a predictive model based on the S U(3)c®S U(3)i <S)U( \ )x gauge group supplemented 
by the A^Z^Z^Z^Zk, discrete group, which successfully describes the SM fermion mass and 
mixing pattern. The small active neutrino masses are generated via inverse seesaw mechanism 
with three very light Majorana neutrinos. The observed charged fermion mass hierarchy and 
quark mixing pattern are originated from the breaking of the Z\®Zp® Z\p discrete group at 
very high scale. The obtained values for the physical observables for both quark and lepton 
sectors are in excellent agreement with the experimental data. The model predicts a vanishing 
leptonic Dirac CP violating phase as well as an effective Majorana neutrino mass parameter of 
neutrinoless double beta decay, with values mpp = 2 and 48 meV for the normal and the inverted 
neutrino mass hierarchies, respectively. 

Keywords: Fermion masses and mixings. Discrete flavor symmetries, 3-3-1 models. Models 
beyond the Standard Model. 


1. Introduction 

Despite the great success of the Standard Model (SM), recently confirmed by the discovery of 
the 126 GeV Higgs boson by LHC experiments MS!,' there are many aspects not yet 
explained such as the origin of the fermion mass and mixing hierarchy as well as the mechanism 
responsible for stabilizing the electroweak scale Si. This discovery of the Higgs scalar field 
allows to consider extensions of the SM with additional scalar fields that can be useful to explain 
the existence of Dark Matter U 

The Standard Model is a theory with many phenomenological achievements. However in the 
Yukawa sector of the SM there are many parameters related with the fermion masses with no clear 
dynamical origin. Because of this reason, it is important to study realistic models that allow to set 
up relations among all these parameters of the Yukawa sector. Discrete flavor symmetries allow 
to establish ansatz that explain the flavor problem, for recent reviews see Refs. mm. These 
discrete flavor symmetries may be crucial in building models of fermion mixing that address the 
flavor problem. Non abelian discrete flavor symmetries arise in string theories due to the discrete 
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features of the fixed points of the orbifolds m . For instance, the discrete D 4 group is originated 
in the S 1 /Z 2 orbifold nn. 

Besides that, another of the greatest mysteries in particle physics is the existence of three 
fermion families at low energies. The quark mixing angles are small whereas the leptonic mixing 
angles are large. Models based on the gauge symmetry S U(3)c®S U(3)l®U(\)x have the feature 


of being vectorlike with three families of fermions and are therefore anomaly free [12,1131, 14 


E3EH- When the electric charge is defined in the linear combination of the S U(3)l generators 
Tt, and 7'x, it is a free parameter, independent of the anomalies (fj). The choice of this parameter 


defines the charge of the exotic particles. Choosing p 


VI 


, the third component of the weak 


lepton triplet is a neutral field v R , which allows to build the Dirac matrix with the usual field v r 
of the weak doublet. If one introduces a sterile neutrino Nr in the model, then it is possible to 
generate light neutrino masses via inverse seesaw mechanism. The 3-3-1 models with /? = --4= 
have the advantange of providing an alternative framework to generate neutrino masses, where 
the neutrino spectrum includes the light active sub-eV scale neutrinos as well as sterile neutrinos 
which could be dark matter candidates if they are light enough or candidates for detection at 
the LHC, if their masses are at the TeV scale. This interesting feature makes the 3-3-1 models 
very interesting, since if the TeV scale sterile neutrinos are found at the LHC, these models can 
be very strong candidates for unraveling the mechanism responsible for electroweak symmetry 
breaking. Furthermore, the 3-3-1 models can provide an explanation for the 750 GeV diphoton 
excess recently reported by ATLAS and CMS [ 17] as we ll a s for the 2 TeV diboson excess found 
by ATLAS [18]. 

Neutrino oscillation experiments fl EH EH 21, EH EH indicate that there are at least two 
massive active neutrinos and at most one massless active neutrino. In the mass eigenstates, it is 
necessary for the solar neutrinos oscillations that 5nr slm = = mi, - m 2 where m 2 ; - m 2 > 0. 

For the atmospheric neutrinos oscillations it is required that 6 m 2 tm = m 2 l = m 2 - m 2 where 
the difference can be possitive (normal hierarchy) or negative (inverted hierarchy). Neutrino 
oscillations do not give information neither on the absolute value of the neutrino mass nor on the 
Majorana or Dirac nature of the neutrino. However there are neutrino mass bounds arising from 
cosmology EH , tritium beta decay EHl and double beta decay 1 26[ 27, 28, EH EH 31,32, 34, 33f] ■ 

The neutrino masses and mixings are known from neutrino oscillations, which depend on 
the squared neutrino mass diferences and not on the absolute value of the neutrino masses. The 
global fits of the available data from the Day a Bay OH. T 2 K EH, MINOS EH, Double CHOOZ 
122] and RENO (EH neutrino oscillation experiments, constrain the neutrino mass squared split¬ 
tings and mixing parameters (35!]. The current neutrino data on neutrino mixing parameters can 
be very well accommodated in the approximated tribimaximal mixing matrix. 


Utbm - 


11 2L 
V3 VI 


1 

V6 

1 

V6 


1 

VI 

1 

VI 


0 

1 

" VI 
1 

V2 


(1) 


which is consistent with two large mixing angles and one very small one of order zero. Specifi¬ 
cally, the mixing angles predicted by the tribimaximal mixing matrix satisfy (sin 2 # 12 ) = 5 , 

(sin 2 023 ) TBM = 5 , an< 4 (sin 2 ®u) TBM = 0- However, the 3-3-1 model is not able to generate the 
tribimaximal matrix structure. Because of this reason, discrete symmetry groups 130. .37. 38.41, 
EH EH EH EE , 44,45] that act on the fermion families are imposed with the aim to generate ansatz 
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that reproduce these matrices. One of the most promising discrete flavor groups is A 4 , since it is 
the smallest symmetry with one three-dimensional and three distinct one-dimensional irreducible 
representations, where the three families of fermions can be accommodated rather naturally. An¬ 
other approach to describe the fermion mass and mixing pattern consists in postulating particular 
mass matrix textures (see Ref [46] for some works considering textures). Besides that, models 
with Multi-Higgs sectors. Grand Unification, Extradimensions and Superstrings as well as with 
horizontal symmetries have been proposed in the literature 
observed pattern of fermion masses and mixings. 


47, 48, |49|, 50] to explain the 


In this paper we propose a version of the S U(3)c x S U(3)l x U (1 )x model with an additional 
flavor symmetry group A 4 ®> Z 3 ®> Z 4 ® ® Z|g and an extended scalar sector needed in order to 

reproduce the specific patterns of mass matrices in the fermion sector that successfully account 
for fermion masses and mixings. The particular role of each additional scalar field and the corre¬ 
sponding particle assignments under the symmetry group of the model are explained in detail in 
Sec. [2] The model we are building with the aforementioned discrete symmetries, preserves the 
content of particles of the minimal 3-3-1 model, but we add additional very heavy scalar fields 
with quantum numbers that allow to build Yukawa terms invariant under the local and discrete 
groups. This generates the predictive and viable textures that explain the 18 physical observables 
in the quark and lepton sectors, i.e., the 9 charged fermion masses, 2 neutrino mass squared split¬ 
tings, 3 lepton mixing parameters, 3 quark mixing angles and 1 CP violating phase of the CKM 
quark mixing matrix. Our model successfully describes the prevailing pattern of the SM fermion 
masses and mixing. 


The content of this paper is organized as follows. In Sec. [2]we outline the proposed model. In 
Sec. [T] we discusss lepton masses and mixings and show our corresponding results. Our results 
for the masses and mixings in the quark sector followed by a numerical analysis are presented in 
Sec. [4] Finally in Sec. 0 we state our conclusions. In |Appendix A| we present a brief description 
of the A 4 group. 


2. The model 

We extend the S U(3)c ® S U(3)i ®> U(l)x group of the minimal 3-3-1 model by adding an extra 
flavor symmetry group A4 ®> Z3 ®> Z4 ®> Z$ <S> Z ig, in such a way that the full symmetry Q experi¬ 
ences a three-step spontaneous breaking, as follows: 

Q = SU(3)c®S U (3 ) l ®>U(l) x ® A 4 ®Z 3 ® Z 4 ®>Z 6 ® Z i6 (2) 

A S U(3) c ®SU(3 ) l ®> U (l) x ® zAs U(3) c ®SU (2) L <g> U (l) y 
-Ast/(3) c ®t/(l) e , 

where the different symmetry breaking scales satisfy the following hierarchy v 7; , v p <s v x <sc A,,,,. 

We define the electric charge in our 3-3-1 model in terms of the SU( 3) generators and the 
identity, as follows: 

Q = T 3 - -^-T & + XI, (3) 

V3 

with I = Diag( 1,1,1), Tt, = \Diag( 1, -1,0) and T% = (^j=)Diag( 1,1,-2). 
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The anomaly cancellation of S U(3)l requires that the two families of quarks be accommodated 
in 3* irreducible representations (irreps). Besides that, the number of 3* irreducible represen¬ 
tations is six, as follows from the quark colors. We accommodate the other family of quarks 
into a 3 irreducible representation. Furthermore, we have six 3 irreps taking into account the 
three families of leptons. Thus, the S U(3)l representations are vector like and free of anomalies. 
Having anomaly free U( 1 )x representations requires that the quantum numbers for the fermion 
families be assigned in such a way that the combination of the U(\)x representations with other 
gauge sectors cancels anomalies. Consequently, to avoid anomalies, the fermions have to be ac¬ 
commodated into the following ( S U(3)c, S U{3 )l, U( 1 )x) left- and right-handed representations: 


' D 1 ’ 2 ' 


ft/ 3 ) 


f v 1 ’ 2 ’ 3 3 

-U 1 ’ 2 

: (3,3*,0), Q 3 L = 

D 3 

: (3,3,1/3), L|’ 2 ’ 3 = 

e l,2,3 

J 1 ’ 2 

L 

J , 

L 

k (v IA3 ) c , 


D 1 / : (3,1,-1/3), Ur : (3,1,2/3), 
U 1 ’ 2 : (3,1,2/3), D\ : (3,1,-i/3), 
Jr : (3,1, -1/3), T r :( 3,1,2/3), 


sr : (1) 1, — 1). /i*:(l,l,-l), r s :(l,l,-l). 

Nr : (1,1,0), Nr : (1,1,0), N>: (1,1,0), 


(1,3,-1/3), 


(4) 


where U' L and D‘ L for i = 1,2,3 are three up- and down-type quark components in the flavor 
basis, while v' L and e' L ( eL,p L , T/J are the neutral and charged lepton families. The right-handed 
fermions are assigned as S U(3)l singlets representations having U( 1 )x quantum numbers equal 
to their electric charges. Furthermore, the fermion spectrum of the model includes as heavy 
fermions: a single flavor quark T with electric charge 2/3, two flavor quarks J 1 ’ 2 with charge 
-1/3, three neutral Majorana leptons (v 1,2,3 )£ and three right-handed Majorana leptons Nr 23 
(see Ref. Jill] for a recent discussion about neutrino masses via double and inverse see-saw 
mechanism for a 3-3-1 model). 

The 3-3-1 models extend the scalar sector of the SM into three 3’s irreps of S U(3)l, where 
one heavy triplet x acquires a vaccuum expectation value (VEV) at the TeV scale, v x , breaking 
the S U(3)l x U(l) x symmetry down to the S U{2)l x f/(l)r electroweak group of the SM and 
then giving masses to the non SM fermions and gauge bosons; and two lighter triplet fields rj and 
p that get VEVs v n and v p , respectively, at the electroweak scale thus generating the mass for the 
fermion and gauge sector of the SM. We enlarge the scalar sector of the minimal 3-3-1 model by 
introducing 14 S U(3)l scalar singlets, namely, £j, S j , ip. A, 0, r and cr (j - 1,2,3). 

The scalars of our model are accommodated into the following [.S' U(3)l, U( \ )x] representa¬ 
tions: 

( „o ^ 

i;i : ( 1 , 0 ), 


£j : (1,0), 


Sj: (1,0), cr ~ (1,0), j= 1,2,3. (5) 


X ■ 


K^ v x + ^x ±l ^J 

P+ 

"^( v p + ^p ± 

Pi 


: (3, -1/3), 


: (3,2/3), 


1 = 


Th 


: (3, -1/3), 


t : (1,0), ^ : (1,0), j= 1,2,3, 

0 : (1,0), A: (1,0), ;= 1,2,3, 
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The scalar fields are grouped into triplet and singlet representions of A 4 . The scalar fields of 
our model have the following assignments under A4 ® Z3 ® Z4 ® Zg ® Z|g: 

77 ~ 1 , 1 , l), p ~ (l,e^, 1 , 1 ,l), x ~ ( 1 . 1 , 1 , 1 , 1 ), 

£ ~ (3,1,1,1,-1), £ ~ (3,1,1,1,e*), S ~(3,e _2 r,i,i, e ¥), cr ~ (l,l, 1, , 

ip ~ (l,l,l,e _ T,l), A~(l,l,-l,e _ T,l), 0~(l',l,i,l,l), r~ (l",l,i, 1,1), (6) 

where the numbers in boldface are dimensions of the A4 irreducible representations. 

The leptons transform under ,4 4 ® Z 3 ® Z 4 ® Z () ® Z\(, as: 


(3,^,1, 1 ,- 1 ), 

e R ~ (l, 1,1, l.e*'), 

1,1,0, 


(l ,l, 1, l,e 7 ) , 

N r ~ (3,e“,l,l,-l). 


(7) 


Note that left handed leptons are unified into a A 4 triplet representation 3, whereas the right 
handed charged leptons are assigned into different A 4 singlets, i.e, 1, V and 1". Furthermore, the 
right handed Majorana neutrinos are unified into a A 4 triplet representation. 

The A 4 ® Z 3 ® Z 4 ® Zg ® Zi 6 assignments for the quark sector are: 


Gi ~ (u,l,l,^), 

1 , 1 , 1 ), 

( 2 i~ (I'M, 1 , 1 , 1 ), 


4 ~ (l,e-x,i,i, e ; r), 

t/| ~ (l\e _2 r,i,i,;' 

), 

.1.1.1). 

D\ ~ (l,e“,l,-l,e¥), 

D* ~(l,e“,l,-l,l 

). Dip IVf, 

1-1.1) 

r s ~ (1",1,1,1,1), 

4 ~ (i',i, 1,1,1) 

4 ~ (I'M, 1,1,1). 



With the above particle content, the following relevant Yukawa terms for the quark and lepton 
sector arise: 
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where the dimensionless couplings y{p , y i '- ) ( i,j = 1,2,3), y (T \ y^\ y^\ h !®, /z®, /z®, /z®, 
/iiw, x, /z 2 ,v and h p are 0( 1) parameters. Here we assumed that all Yukawa couplings are real, 
excepting y^\ y^ and /z® which are assumed to be complex. 

Although the flavor discrete groups in Eq. (0 look rather sofisticated, each discrete group 
factor plays its own role in generating predictive fermion textures that successfully account for 
the low energy fermion flavor data. To describe the pattern of fermion masses and mixing angles, 
one needs to postulate particular Yukawa textures. As we will see in the next sections, the 
predictive textures for the lepton and quark sectors will give rise to the experimentally observed 
deviation of the tribimaximal mixing pattern and to quark mixing emerging only from the down 
type quark sector, respectively. A candidate for generating specific Yukawa textures is the A 4 
flavor symmetry, which needs to be supplemented by the Z 3 ®Z 4 ®Zg®Zi 6 discrete group. As we 
will see in the next sections, this predictive setup can successfully account for fermion masses and 
mixings. The inclusion of the A 4 discrete group reduces the number of parameters in the Yukawa 
and scalar sector of the S U(3)c 8 > S 1/(3)/, ® U( l)x model making it more predictive. We choose 
A 4 since it is the smallest discrete group with a three-dimensional irreducible representation and 
3 distinct one-dimensional irreducible representations, which allows to naturally accommodate 
the three fermion families. We unify the left-handed leptons in the A 4 triplet representation 
and the right-handed leptons are assigned to A 4 singlets. Regarding the quark sector, we assign 
quarks into A 4 singlet representations. In what follows we describe the role of each discrete 
cyclic group factor introduced in our model. The Z 3 symmetry separates the A 4 scalar triplets 
participating in the Yukawa interactions for charged leptons from those ones participating in the 
neutrino Yukawa interactions. Besides that, the Z 3 symmetry avoids mixings between SM quarks 
and exotic quarks since the right handed exotic quarks are neutral under this symmetry whereas 
the right handed SM quarks have non trivial Z 3 charges. Thus the Z 3 symmetry decouples the 
SM quarks from the exotic quarks resulting in a reduction of quark sector model parameters. 
Furthermore, the Z 4 symmetry is also important for reducing the number of quark sector model 
parameters, since due to this symmetry, the S U(3)i scalar singlets A 4 nontrivial singlets only 
appear in the down type quark Yukawa terms. Consequently this Z 4 symmetry together with 
the A 4 assignments for quarks described in Eq. 0, results in a diagonal up type quark mass 
matrix, thus giving rise to a quark mixing only emerging from the down type quark sector. The 
Z(, symmetry is crucial for explaining the hierarchy between the SM down and SM up type 
quarks without tuning the SM down type quark Yukawa couplings, since it is the smallest cyclic 
symmetry that allows p in the Yukawa term that generates the bottom quark mass, which is T 3 -jj= 
(A = 0.225 is one of the Wolfenstein parameters) times a 0(1) parameter. The Z 16 symmetry 
gives rise to the observed hierarchy among charged fermion masses and quark mixing angles. 
It is worth mentioning that the properties of the Zy groups imply that the Z 16 symmetry is the 

_ 1 g 

smallest cyclic symmetry that allows to build the Yukawa term Q L p*U R ^j of dimension twelve 

from a p insertion on the (jlfA U r operator, crucial to get the required A H suppression (where 
A = 0.225 is one of the Wolfenstein parameters) needed to naturally explain the smallness of the 
up quark mass. Regarding the charged lepton sector, let us note that the five dimensional Yukawa 
operators j ( ' L R pe R , ^ ( Lipp R and ^ ( L R p^j t r are A 4 invariant but do not preserve the 
Z 16 symmetry, as follows from the charges assignments given by Eqs. © and 0 . 

In what follows we comment about the possible VEVs patterns for the A 4 scalar triplets £ 
and S. Here we assume a hierarchy between the VEVs of the A 4 scalar triplets A zT and S , i.e., 
v s « \’i « vy , which implies that the mixing angles of these scalar triplets are very small 
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since they are suppressed by the ratios of their VEVs, which is a consequence of the method of 
recursive expansion proposed in Ref. (52]. Consequently, we can neglect the mixing between 
the A\ scalar triplets 1 and S , and treat their corresponding scalar potentials independently. 
The relevant terms determining the VEV directions of any A 4 acalar triplet are: 

V (Z) = -a4 (ZZ*)i + ATS,! (EE*)! (ZZ*)j + Kta (SZ), (Z*Z*)! + * 1,3 (ZZ*) r (ZZ*)i„ 

+*x ,4 [(ZZ) r (Z*Z*) r , + h.c] + ^,5 [(ZZ),„ (Z*Z*) F + h.c] 

+KZ ,6 (ZZ*) 3s (ZZ*) 3s + K'ej (ZZ) 3s (Z*Z*) 3s . (11) 

where Z = g, g, S. 

The part of the scalar potential for each A 4 scalar triplet has 8 free parameters: 1 bilinear and 7 
quartic couplings. The minimization conditions of the scalar potential for a A 4 triplet yield the 
following relations: 

-2vx,/4 + 4 Kr A \ % (vi, + vl 2 + v| 3 ) + 2 *s, 3 vs 1 (2v^ - v^ - v§ 3 ) 

+4KI.2VX, [vi, + v| 2 cos ( 201 , - 20s, ) + v| 3 cos (20s, - 20s,)] + 8/^sjVs, (v| 2 + v| 3 ) 
+4 («s,4 + * 5 , 5 ) Vs, [2v| 3 - v| 2 cos (20s, - 20 s 2 ) - v| 3 cos (20s, - 20s 3 )] 

+4a-s, 6 vs, [v|, {1 + cos (20s, - 20v,)) + v| 3 {1 + cos (20s, - 20v 3 )}] 

0, 

-2vs 2 //s + (vl, + vi, + vi 3 ) + 2 *s, 3 vs 2 (2vi 2 - vi, - vf 3 ) 

+4a-s, 2 vs, [vi 2 + vi i cos (20s, - 20 E ,) + vi 3 cos (20v, - 20s,)] + 8ksjVs 2 (vi, + vi 3 ) 
+4 (/ts,4 + * 1 , 5 ) vs 2 [2t'i, - vi, cos (20v, - 20s,) - v| 3 cos (20s, - 20s 3 )] 

+4a-s, 6 vs, [v|, {1 + cos (20v, - 20s,)} + i'i 3 {1 + cos (20s, - 20v 3 )}] 

0, 

-2vz 3 Ati + 4/ts,tVs 3 (vi, + vi, + vi 3 ) + 2k X3 vv 3 (2vi 3 - vi, - vfj 

+4a-s, 2 vs 3 [vi, + vi, cos (20s, - 20 x ,) + i'i, cos (20s 3 - 20s,)] + 8 ksjVs 3 (vi, + vi,) 
+4 (/ts,4 + * 1 , 5 ) vs 3 [2vi 3 - vi, cos (20s, - 20s,) - vi, cos (20s 3 - 20s,)] 

+4a"s,6Vs 3 [vi, {1 + cos (20s, - 20s,)} + vi, {1 + cos (20s 3 - 20s,)}] 

0. (12) 

where (Z) = (vs,e' Sx i, \%e ,BL 2 , vs^" 9 ^). Here in order to simplify the analysis, we restrict to the 
simplest case of zero phases in the VEV patterns of the A 4 triplet scalars, i.e., 0s, = 0s, = 0s 3 = 0. 
Then, from the scalar potential minimization equations given by Eq. (fl2T> . the following relations 
are obtained: 


d<V(Z)) 

d\% 


d<V(Z)) 

dvs. 


d<V(Z)) 

di's 3 


[3/Cs,3 - 4 (k-£,6 + + 6 (ks.4 + ^, 5)] (vi, - l'i,) = 

[3ats,3 - 4 (ks,6 + *1,7) + 6 (*s,4 + *£,5)] (vi, - Vv 3 ) = 

[3/CS.3 - 4 (/Cs.6 + Ky,i) + 6 (/Cs,4 + ^, 5 )] (vi 2 - vi 3 ) = 

7 


0 , 

0 , 

0 . 


( 13 ) 





From the relations given by Eq. (IT3l > and setting 3 = | + ^, 7 ) - 2 (^,4 + ^, 5 ). we obtain 

the following VEV pattern: 


<f> = -^=(1,1,1), 
V3 


(0 = —p ( 1 , 0 , l), 
V 2 


(S) = —= (1,1,-1) • 

V3 


(14) 


) or / 14 , which has been 
s. SIM)- The VEV 


In the case of this is a vacuum configuration preserving a Z 3 subgroup of ,4 4 
extensively studied in many A 4 flavor models (for recent reviews see Refs, 
pattern for the A 4 triplet scalar C, is similar to the one previously studied in an A 4 and Tj flavor 
S U( 5) GUT models [38’, |44] and in a 6 HDM with A 4 flavor symmetry 137]. As we will see in 
the next section, the VEV patterns for the A 4 triplets £, £ and S given in Eq. IT41) are crucial 
to get a predictive model that successfully reproduces the experimental values of the physical 
observables in the lepton sector. 


Furthermore we assume that these S U(3)l scalar singlets get VEVs at a scale K int much larger 
than v x (which is of the order of the TeV scale), with the exception of S ; (j =1,2,3), which get 
VEVs much smaller than the electroweak symmetry breaking scale v = 246 GeV. The VEVs of 
the gj ( j =1,2,3), ip. A, (f>, t and cr scalar singlets break the S U(3)c ® S U(3)l ® U(l)x ® A 4 ® 
Z3 ® Z 4 ® Zg ® Zi6 symmetry down to S U(3)c ® S U(3)l ® U(l)x ® Z3 at the scale A 


From the expressions given above, and using the vacuum configuration for the A 4 scalar triplets 
given in Eq. (fl4t . we find the following relations: 


£ = 

2 

| [3 (/c f ,| + K(, 2 ) + 4 (/Cf,6 + K f , 7 )] v]. 


$ = 

2 

- [3 (k (a + Kp t 2 ) + 4 (k l6 + /C f , 7 )] v\. 


Vs = 

2 , , 

2 P fc.i + «s, 2 ) + 4(/c 5 ,6 + *s, 7 )] v s . 

(15) 


These results show that the VEV directions for the three A 4 triplets, i.e., C and S scalars in Eq. 
<ns, are consistent with a global minimum of the scalar potential (ITTI) of our model for a large 
region of parameter space. 

Besides that, as the hierarchy among charged fermion masses and quark mixing angles emerges 
from the breaking of the Z 4 ® Zg ® Zi6 discrete group, we set the VEVs of the S U(3)l singlet 
scalar fields 4 ; , <p. A, <p, r and cr, with respect to the Wolfenstein parameter A = 0.225 and the 
model cutoff A, as follows: 


~ V T ~ ~ V A ~ Vf ~ Vo- ~ K in1 = TA. (16) 

Furthermore, we assume that the A 4 scalar triplets participating in the neutrino Yukawa inter¬ 
actions have VEVs much smaller than the electroweak symmetry breaking scale. Besides that, 
as previously mentioned, we assume a hierarchy among the VEVs of the two A 4 scalar triplets 
participating in the neutrino Yukawa terms. Consequently, as we will see in the next section, 
the Majorana neutrinos acquire very small masses and thus an inverse seesaw mechanism for 
the generation of light active neutrino masses, takes place. Therefore, we have the following 
hierarchy among the VEVs of the scalar fields in our model: 


Vs « Vf « V p ~ v n ~ V « v x « K m ,. 


(17) 
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In what follows, we briefly comment about the low energy scalar sector of our model. The 
renormalizable low energy scalar potential of the model is given by: 

Vh = H 2 X (X*X) + + Pp(p f p) + f {xiTljPk^ Jk + He •) + M(x f X)(X f X) 

+/l 2 (p t p)(P t P) + + M(X } X)(P*P) + 'hix^xW’l) 

+A 6 (p f p)(ifTj) + A^TfArfx) + ^s( X f P)(P f X) + ^9^rj)(jfp). (18) 


After the symmetry breaking takes place, it is found that the scalar mass eigenstates are related 
with the weak scalar states by: [14, 15]: 



= R. 


‘( J ‘T 





with 


/ cos a T (/3 T ) sin a T (J3 T ) \ 

\ - sin a T (/3 T ) cos a T (/3 T ) )’ 



where tan [3 T = v^/vp, and tan 20 : 7 - = M\/(M 2 - M 3 ) with: 


(19) 

( 20 ) 


( 21 ) 


M, 

= 4/l6V,v' p + 2 V2 fv x . 


m 2 

- 44 2 V 2 - V2/v a tan/l r , 


m 3 

= 4/t 3 v^ - ^/lfv x / tan/l-p. 

(22) 


It is noteworthy to mention that the our model has the following scalar states at low energies: 4 
massive charged Higgs (Hf, H 2 ), one CP-odd Higgs (A®), 3 neutral CP-even Higgs (h°, HH { \) 

and 2 neutral Higgs bosons. We identify the scalar h° with the SM-like 126 GeV Higgs 

boson discovered at the LHC. Let us note that the neutral Goldstone bosons Gj, G®, G 2 , G 2 corre¬ 
spond to the longitudinal components of the Z, Z', K {] and K gauge bosons, respectively. Besides 
that, the charged Goldstone bosons Gj and Gj are associated to the longitudinal components of 
the W ± and K 1 gauge bosons, respectively [ 12|, 15]. 
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3. Lepton masses and mixings 


From Eq. (fTOl) and taking into account that the VEV pattern of the A 4 triplet, S U(3)l singlet 
scalar field £ satisfies Eq. (IT4l) with the nonvanishing VEVs of the S U(3)i singlet scalars £ and 
cr, set to be equal to A A (being A the cutoff of our model) as indicated by Eq. we find that 
the charged lepton mass matrix is given by: 


(), R/ L = - 7 = 

( 

1 

1 

1 


' 1 

0 

0 ' 


1 

CJ 

co 2 

2 m 

, lj — e 3 , Pi — 

0 

1 

0 

V H / -/T 

V 

1 

OJ 2 

V , 


l 0 

0 

e m , 


being a the complex phase of hi 1 /, and the charged lepton masses are: 


(23) 


m e — a'lA s -, 

V 2 


nia = at? A 5 ——, 

V 2 


m T = a? A 3 

V2 


(24) 


where A — 0.225 is one of the Wolfenstein parameters, v = 246 GeV the scale of electroweak 
symmetry breaking and a l> (i = 1,2,3) are 0(1) parameters. Let us note that the charged lepton 
masses are linked with the scale of electroweak symmetry breaking through their power depen¬ 
dence on the Wolfenstein parameter A = 0.225, with 0(1) coefficients. Furthermore, from the 
lepton Yukawa terms given in Eq. dTOt it is easy to see that our model does not feature flavor 
changing leptonic neutral Higgs decays. Consequently, our model cannot explain the recently 
reported anomaly in the h —» fir decay, implying that a measurement of its branching fraction 
will be decisive for its exclusion. 

Regarding the neutrino sector, we can write the neutrino mass terms as: 



( VL } 

( yf V R Nr )M v 

V R 

l < J 


+ H.c, 


(25) 


where the neutrino mass matrix is constrained from the A 4 flavor symmetry and has the following 
form: 

03x3 M d 03x3 N 

M t d 0 3 x 3 M x , (26) 

0 3 x3 Ml Mr 


My = 


and the submatrices are given by: 


M d 


M x 


hpVpV^Va- 

2A 1 


0 

-1 

0 


1 0 ' 

0 -1 

1 0 J 


1 0 0 f 
0 1 0 
^2 0 0 1 



( 1 '’l +xv l 

h\N A 2 


U V S V <r ) 

,l2N HR 

Mr = 

U v sV(T 

,l2N V3A 

hm^-Va- 

U V <r 

h2N VfX 


U v s v <y 

< - N 1E 

U v SV<r 

,l2N lfX 

"IN A 2 V(T / 


(27) 


As previously mentioned, we assume in our model that the S U (3)/, scalar singlet, A 4 triplet S 
interacting with the right handed Majorana neutrinos gets a very small vacuum expectation value, 
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much smaller than the electroweak symmetry breaking scale, which results in very small masses 
for these Majorana neutrinos. Consequently, this setup can generate small active neutrino masses 
through an inverse seesaw mechanism. 

As shown in detail in Ref. mb, the full rotation matrix is approximately given by: 


U = 


Ry 

V2 

(. B\-B \) 

V2 


B 3 U x 

fiz Sin 

V5 u x 
tlzS }tj 
V2 V X 


B 2 Ur 
(l+S) 

V2 
(1-S) 


V2 


Ur 

Ur 


where 


.S' = - 


2 V 2 hfv x 


Mr, 


and the physical neutrino mass matrices are: 


B 2 — B 3 — 




-M 


D’ 


M R M x l M T D , 

m[ 2) = -m t x +\m r , m (3) = m t x + 1 -m r . 


(28) 


(29) 


(30) 

( 31 ) 


where M\ X) corresponds to the active neutrino mass matrix whereas m[ 2> and m[ 3) are the exotic 
Dirac neutrino mass matrices. Note that the physical eigenstates include three active neutrinos 
and six exotic neutrinos. The exotic neutrinos are pseudo-Dirac, with masses ~ ±M^ and a small 
splitting Mr. Furthermore, R v , Ur and U x are the rotation matrices which diagonalize M\}\ m[ 2) 


and M v , respectively [51]. 


From Eq. (l30l > it follows that the light active neutrino mass matrix is given by: 


M (2) 


A 


h 2 pV 2 p v 2 Vo- 


/ V 2 

hm± 

0 

+ 

hm'f 



' A 

0 

A ' 

0 

I, v x , 2/ijw 

h lN a + ~jfV S 

0 

0 


= 

0 

B 

0 

. (iuv-f 


7 V X 

> 


l A 

0 

^4 J 


2^v?A4 


B = — 


n p V p V l v o 


ihf'vlrs 


h in~t + 
A 


2 h 2 N 

V3 


,, 2 \ 


+ /tnv-f 
A 


(32) 


The neutrino mass matrix given in Eq. (l32l > only depends on two effective parameters: A and 
B. These effective parameters include the dependence on the various model parameters. It is 
noteworthy that A and B are suppressed by inverse powers of the high energy cutoff A of our 
model. 

The light active neutrino mass matrix M[ X) is diagonalized by a unitary rotation matrix R v , 
according to: 


RyM ( v l) R v = 


777 1 

0 

0 ' 


cos 8 

0 

sin6> ' 

0 

777 2 

0 

, with R v = 

0 

1 

0 

o 

0 

777 3 , 


_ - sin 8 

0 

COS0 , 


71 

* = ± 4’ 


(33) 


where the upper sign corresponds to normal (8 = +7r/4) and the lower one to inverted (8 = — tt/ 4) 
hierarchy, respectively. The light active neutrino masses for the normal (NH) and inverted (IH) 
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mass hierarchies are given by: 


NH 

: 0 = 

11 

+ 

-M 

m Vl = 0 , 

m Vl = B, 

m V3 = 2A, 

(34) 

IH 

: 6 ~- 

I't- 

1 

m Vl = 2 A, 

m Vl = B, 

»z v , = 0 . 

(35) 


We also find that the PMNS leptonic mixing matrix is given by: 


vz/ 


cosg 

VI 

cos 6 

VI 


sin 9 ia 

V3 g 


sin 9 ia+- 
V3 g 


1 

V3 


cos v ia 
V3 g 


sing 

V3 


V3 


e 3 


cos 6 _ sin# 1 0 = 1 t- 

V3 V3 g V3 C 


cos 8 „ia+ Mr- , sing 

V3 g + VI 


cos 6 „ia- =fr- , sing 

VI e + VI 


(36) 


It is worth commenting that the Pontecorvo-Maki-Nakagawa-Sakata(PMNS) mixing matrix de¬ 
pends only on the parameter a, while the neutrino mass squared splittings are controlled by 
parameters A and B. 

The standard parametrization of the leptonic mixing matrix implies that the lepton mixing angles 
satisfy (<§]: 


--, Sin 2 013 = \U e i\ 2 = i(l + cos a), (37) 

cos a 3 

V3 sin a ' 

- + -^ ■ 

cos a — 2) 

The resulting PMNS matrix (l36t reduces to the trimaximal mixing matrix O in the limit a = n, 
for the inverted and normal hierarchies of the neutrino mass spectrum. Let us note that the 
lepton mixing angles are controlled by a single parameter (a), whereas the neutrino mass squared 
splittings only depend on the parameters A and B. 

The Jarlskog invariant and the CP violating phase are 0: 


sin" 0 i 2 = 


\Ue2? 


l-\U e3 \ 2 2- 

2 


sin 2 #23 = 


\M 


i 


1-1 U e3 \ Z 


J = hn(u el U ia U‘ a U; i ) 


-cos 20, 

6 V3 


sin <5 = 


_ 8 J _ 

cos 0[3 sin 20n sin 2023 sin 2013 


■ (38) 


Taking into account that 0 = +1, our model predicts J = 0 and <) = 0, which results in a vanishing 
leptonic Dirac CP violating phase. 

In what follows we adjust the three free effective parameters a, A and B of the lepton sector of 
our model to reproduce the experimental values of the five physical observables in the neutrino 
sector, i.e., three leptonic mixing parameters and two neutrino mass squared splittings, reported 
in Tables [DEI for the normal (NH) and inverted (IH) hierarchies of the neutrino mass spectrum, 
respectively. We fit only a to adjust the values of the leptonic mixing parameters sin 2 0, ; , whereas 
A and B for the normal (NH) and inverted (IH) mass hierarchies are given by: 


NH : m Yl = 0, 


= B = 



« 9meV, 


m V3 — 2 A - 



« 51meV; 


(39) 
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IH : m Vl = B = Am 2 } + Am 2 3 « 50meV, m Vl — 2A — ^Ja m 2 n ~ 49meV, m r3 = 0, (40) 

as resulting from Eqs. dS),® and the definition Am 2 . = m 2 - m 2 . The best fit values of Am 2 . 

I_- __ l J t J IJ 

have been taken from Tables|T|and[2]for the normal and inverted mass hierarchies, respectively. 

We vary the model parameter a in Eq. (l37l > to fit the leptonic mixing parameters sin 2 0, ; - to the 
experimental values reported in Tables |T] [I] We obtain the following best fit result: 


NH : 

: a = -0.887T, 

sin 2 0 i 2 ~ 0.34, 

sin 2 023 w 0.61, sin 2 0i3 » 0.0232; 

(41) 

IH 

: a = 0.12 7r, 

sin 2 0 i 2 ~ 0.34, 

sin 2 023 ~ 0.61, sin 2 0 |3 =» 0.0238. 

(42) 


From the comparison of Eqs. (l42l) . (ITO with Tables [HO it follows that sin 2 $13 and sin 2 023 
are in excellent agreement with the experimental data, for both normal and inverted mass hierar¬ 
chies, whereas sin 2 f )\2 is deviated 2<x away from its best fit values. This shows that the physical 
observables in the lepton sector obtained in our model are consistent with the experimental data. 
Furthermore, as previously mentioned, our model predicts a vanishing leptonic Dirac CP violat¬ 
ing phase. 


Parameter 

Amj, (10 5 eV 2 ) 

Amjj(10 3 eV 2 ) 

(sin 2 0 12 ) 

\ /exp 

(sin 2 023 ) 

\ /exp 

(sin 2 0i 3 ) 

V /exp 

Best fit 

7.60 

2.48 

0.323 

0.567 

0.0234 

lcr range 

7.42 - 7.79 

2.41-2.53 

0.307-0.339 

0.439 - 0.599 

0.0214-0.0254 

2<t range 

7.26-7.99 

2.35-2.59 

0.292-0.357 

0.413-0.623 

0.0195-0.0274 

3<t range 

7.11-8.11 

2.30-2.65 

0.278-0.375 

0.392 - 0.643 

0.0183-0.0297 


Table 1: Range for experimental values of neutrino mass squared splittings and leptonic mixing parameters, taken from 
Ref. f35ll . for the case of normal hierarchy. 


Parameter 

Am 2 j(10 -5 eV 2 ) 

Am 2 3 (10“ 3 eV 2 ) 

(sin 2 0i 2 ) 

\ /exp 

(sin 2 023 ) 

\ /exp 

(sin 2 0o) 

\ /exp 

Best fit 

7.60 

2.38 

0.323 

0.573 

0.0240 

lcr range 

7.42 - 7.79 

2.32-2.43 

0.307-0.339 

0.530-0.598 

0.0221 -0.0259 

2cr range 

7.26-7.99 

2.26-2.48 

0.292-0.357 

0.432-0.621 

0.0202 - 0.0278 

3cr range 

7.11-8.11 

2.20-2.54 

0.278-0.375 

0.403 - 0.640 

0.0183-0.0297 


Table 2: Range for experimental values of neutrino mass squared splittings and leptonic mixing parameters, taken from 
Ref. E3I . for the case of inverted hierarchy. 

In the following we proceed to determine the effective Majorana neutrino mass parameter, which 
is proportional to the amplitude of neutrinoless double beta (0 v/3/3) decay. This effective Majorana 
neutrino mass parameter has the form: 


mpp = 


Yj U ek m n 


(43) 


where U 2 . and m Vk are the PMNS mixing matrix elements and the Majorana neutrino masses, 
respectively. 
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Using Eqs. (l36l >. (l39t . < 14()b and (l43t . it follows that the effective Majorana neutrino mass param¬ 
eter, for both Normal and Inverted hierarchies, acquires the following values: 


2 meV for NH 
"W ~ 1 47 meV f or IH 


(44) 


Our results for the effective Majorana neutrino mass parameter given above, are beyond the 
reach of the present and forthcoming Ov/3/3 decay experiments. The EXO-200 experiment [26j] 
sets the current best upper limit on the effective neutrino mass parameter equal to mpp < 160 
meV, correspoding to > 1.6 x 10 25 yr at 90% C.L. This bound will be improved 


1 / z i—| p"L 

within the not too distant future. The GERDA “phase-II”experiment 127L128 1 is expected to reach 
T^I^C 6 Ge) > 2 x 10 26 yr, which corresponds to mpp <100 meV. A bolometric CUORE exper¬ 
iment, using 130 Te l29tl . is currently under construction. This experiment features an estimated 
sensitivity of about ( VM] 'Vt) ~ 10 26 yr, corresponding to an effective Majorana neutrino 
mass parameter mpp < 50 meV. Besides that, there are proposals for ton-scale next-to-next gen¬ 
eration Ov/3/3 experiments using 136 Xe [30| 221 and 76 Ge 127J, 320, which claim sensitivities over 

rpOvpp 
1/2 


10“' yr, corresponding to an effective Majorana neutrino mass parameter mpp ~ 12 - 30 


meV. For a recent review, see for example Ref. 83411 . Consequently, the Eq. ( 144b indicates 
that our model predicts at t ^ le l eve l of sensitivities of the next generation or next-to-next 
generation Ov/3/3 experiments. 


4. Quark masses and mixing. 


From the quark Yukawa terms of Eq. (0 and the relation given by Eq. (IT6l >. we find that the SM 
quarks do not mix with the exotic quarks and that the SM quark mass matrices are: 


Mu = 


a\ U) A & 


0 

0 


0 

4°* 

0 


0 

0 

J.V) 


v 

1? 


-=, M d = 


(D) ,7 

11 

(£»,6 


a/’A! 
J D) 36 . 
-*31 


« 12^ 6 

«32^ 5 


a^ > A (, e~ is ‘ l 

a ®/! 5 

a^A 3 


\ ( JO) 37 JO) ,6 JO) \ 

V 

v!’ 

(45) 

where A = 0.225 is one of the Wolfenstein parameters, v = 246 GeV the scale of electroweak 
symmetry breaking and aj. 77 ' 0 * 0 , j = 1 , 2 ,3) are 0 ( 1 j parameters. 

Moreover, we find that the exotic quark masses are: 


mj — y 


_ JT) v r 


V2 


m J' = >’, 


_ JJ) _ y ) 


,(J) 


V2 


m 


rrir. 


(J) 

( 7 ) v x y_i 

nip - y\ — = — —m T . (46) 

V2 y (T) 


Since the charged fermion mass and quark mixing pattern emerges from the breaking of the 
Z\ ® Zp ® Zip discrete group and in order to simplify the analysis, the following scenario is 
considered: 

arg («i 3 ) ) = arg , a ^ = a'f, i, j = 1,2,3. (47) 

Besides that, to show that the quark textures given above can fit the experimental data, and in 
order to simplify the analysis, we adopt a benchmark where we set = af ) = 1 and t/!/,’ = 
a 7/7) , as suggested by naturalness arguments and by the relation m c ~ nib , respectively. Then, we 
proceed to fit the parameters a ^ a ( £\ a^\ a { ®\ a ( ®\ a ^ and the phase 6 q , to reproduce the 10 
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physical observables of the quark sector, i.e., the six quark masses, the three mixing angles and 
the CP violating phase. The obtained values for the quark masses, the three quark mixing angles 
and the CP violating phase 6 in Table[3]correspond to the best fit values: 


(D ) 
n 

u n 

=* 1 . 11 , 

(D ) 
n 

u 22 

=* 0 . 59 , 

a {D) = 
U 12 

* 0.54, 



(D ) 

n 

= 0 . 43 , 

(D ) 
n 

u 23 

=* 1.13, 

a (D) = 
“33 

* 1.42, 

6 q == 66°. 

(48) 


Observable 

Model value 

Experimental value 

m u (MeV) 

1.14 

1 -0.45 

m c (MeV ) 

635 

635 ± 86 

m,(GeV ) 

173.9 

172.1 + 0.6 + 0.9 

m c ){MeV) 

2.9 

7 q+u.s 

■^—0.4 

m s (MeV) 

57.7 

C 7 7 +I 6.8 

J 7 -15.7 

nib(GeV) 

2.82 

9 09 + 0 . 1)9 

Z " OZ< -0.04 

sin 0 i 2 

0.225 

0.225 

sin 023 

0.0412 

0.0412 

sin 0 i 3 

0.00352 

0.00351 

6 

66 ° 

68 ° 


Table 3: Model and experimental values of the quark masses and CKM parameters. 


The obtained quark masses, quark mixing angles and CP violating phase exhibit an excellent 
agreement with the experimental data. Let us note, that despite the aforementioned simplifying 
assumptions that allow us to eliminate some of the free parameters, a good fit with the low 
energy quark flavor data is obtained, showing that our model is indeed capable of a very good 
fit to the experimental data of the physical observables for the quark sector. The obtained and 
experimental values for the physical observables of the quark sector are reported in Table [3] We 
use the experimental values of the quark masses at the Mz scale, from Ref. [[53] (which are 
similar to those in I5410 . whereas the experimental values of the CKM parameters are taken from 
Ref. i. 

In what follows we briefly comment about the phenomenological implications of our model in 
the flavor changing processes involving quarks. As previously mentioned, the different Z 3 charge 
assignments for SM and exotic right handed quark fields imply the absence of mixing between 
them. Due to the absence of mixings between SM and exotic quarks, the exotic T , /' and J 2 
quarks do not exhibit flavor changing neutral decays into SM quarks and gauge bosons, SM light 
126 GeV Higgs boson and SM quarks. Thus, assuming that the H® and H 2 neutral Higgs bosons 
are heavier than the exotic T, /' and / 2 quarks, it follows that the flavor changing neutral exotic 
quark decays are absent in our model. Consequently these exotic quarks can be searched at the 
LHC via their flavor changing charged decays into SM quarks and gauge bosons, specifically in 
their dominant decay modes T —> bW and J 1,2 —> tW. These exotic quarks can be produced 
at the LHC via Drell-Yan proccesses mediated by charged gauge bosons, where the final states 
will include the exotic T quark with a SM down type quark as well as any of the exotic J or J 2 
quarks with a SM up type quark. Furthermore, from the quark Yukawa terms, one can easily see 
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that the our model predicts the absence of flavor changing top quark decays t —> he and t —> hit 
at tree level. The flavor changing top quark decays t —> he and t —> hit only arise at one loop 
level and will involve virtual charged gauge bosons and exotic quarks running in the loops. Thus, 
a measurement of the branching fraction for the t —> he and t —> hit decays at the LHC will be 
crucial for confirming or ruling out our model. It would be interesting to perform a detailed study 
of the exotic quark production at the LHC, the exotic quark decay modes and the flavor changing 
top quark decays. This is beyond the scope of this work and is left for future studies. 


5. Conclusions 

We constructed a predictive S U(3)c ® S U(3)i ®> U(l)x model with /? = —^=, based on the A 4 
flavor symmetry supplemented by the Z 3 ®> Z 4 ®> Zp ® Z 16 discrete group. Our model success¬ 
fully accounts for the observed fermion masses and mixing angles. The obtained values for the 
physical observables in both quark and lepton sectors exhibit an excellent agreement with the ex¬ 
perimental data. The A 4 , Z 4 and Z 3 symmetries allow to reduce the number of parameters in the 
Yukawa terms, increasing the predictivity power of the model. The breaking of the Z 4 ®Zp ® Z ;6 
discrete group at high energy, gives rise to the observed charged fermion mass pattern and quark 
mixing hierarchy. In our model the Majorana neutrinos acquire very small masses, much smaller 
than the Dirac neutrino masses, thus giving rise to an inverse seesaw mechanism for the gen¬ 
eration of the light active neutrino masses. In this scenario, the spectrum of neutrinos includes 
very light active neutrinos and TeV scale pseudo Dirac nearly degenerate sterile neutrinos. Our 
model predicts a vanishing leptonic Dirac CP violating phase as well as an effective Majorana 
neutrino mass, relevant for neutrinoless double beta decay, with values mpp — 2 and 48 meV, for 
the normal and the inverted hierarchies, respectively. For the inverted hierarchy neutrino mass 
spectrum, our obtained value of 48 meV for the effective Majorana neutrino mass is within the 
declared reach of the next generation bolometric CUORE experiment (29j or, more realistically, 
of the next-to-next generation tone-scale Ov/i/i-decay experiments. Under the assumption that 
the exotic T, J l and / 2 quarks are lighter than the /-/!,’ and H 2 neutral Higgs bosons, our model 
predicts the absence of the flavor changing neutral exotic quark decays, which implies that they 
can be searched at the LHC via their dominant flavor changing charged decay modes T —> bW 
and 7 1,2 -4 tW. Furthermore, our model predicts the absence of flavor changing neutral top 
quark decays at tree level, implying that they occur at one loop level. Possible directions for 
future work along these lines would be to study the constraints on the heavy charged gauge bo¬ 
son masses in our model arising from the upper bound on the branching fraction for the flavor 
changing top quark decays, the oblique parameters, the Zbb vertex and the Higgs diphoton signal 
strength. The heavy exotic quark decays and their production at the LHC may be useful to study. 
All these studies require carefull investigations that we left outside the scope of this work. 
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Appendix A. The product rules for A 4 

The A 4 group has one three-dimensional 3 and three distinct one-dimensional l, V and 1" irre¬ 
ducible representations, satisfying the following product rules: 

3®3 = 3 J ©3 fl ffilffil'©l", (A.l) 

1®1 = 1, 1'®1" = 1, 1'®1' = 1", 1 " ® 1 " = V, 

Considering (xi,yi,zi) and {x 2 ,y 2 ,z 2 ) as the basis vectors for two Azi-triplets 3, the following 
relations are fullfilled: 


(3® 3)! = xiyi + x 2 y 2 + x 3 y 3 , 

(3 ® 3) 3j = (. x 2 y 3 + x 3 y 2 , x 3 yi + x\y 3 ,x\y 2 + x 2 y\), 
(3 ® 3) 3o = (x 2 y 3 - x 3 y 2 , x 3 y 1 - x\y 3 , x\y 2 - x 2 y 1 ), 


(A.2) 

(3 ® 3) r = xyyi + iox 2 y 2 + oj 2 x 3 y 3 , 
(3 ® 3V = xiyi + arx 2 y 2 + a>x 3 y 3 . 


where oj = e l > . The representation 1 is trivial, while the non-trivial V and 1" are complex 
conjugate to each other. Some reviews of discrete symmetries in particle physics are found in 
Refs 


u gat e to eacn ot 

.HBdUE!. 
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